
% Estimando 'e' con Maxwell-Boltzman.

% Cantidad de celdas:
C = 100;
% Cantidad de bolas:
B = 50;

% Cantidad de simulaciones:
N = 100;

while 1
    % Vector con las respuestas a ¿La primera celda está vacía?
    vacia = [];

    for k = 1:N
        % Se distribuyen al azar las bolas en las celdas:
        estados = ceil(rand(B, 1)*C);

        % ¿La celda 1 está vacía?
        vacia(end+1) = (sum(estados == 1) == 0);
        %        _____/                  \____
        %       | ¿Cuántas hay en la celda 1? |
        %       -------------------------------
    end

    % Estimación de 'e':
    mean(vacia)^(-C/B)
    pause(0.25)
end
